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Qj '. Abstract 



K.^ The variation of current / with voltage V for poly(phenylene vinylene) and 

Qv^ ' other polymer light-emitting diodes has been attributed to carriers tunneling 

o' 

CN ' into broad conduction and valence bands. In actuality the electrons and holes 



tunnel into polaron levels and transport is by hopping among these levels. We 



C^ ' show that for small injection the I-V characteristic is determined mainly by 

a" 

the image force, for large injection by space charge effects, but in both cases 



the strong variation of mobility with field due to disorder plays an important 
role. 



k> ; PACS: 73.40.Ns,72.20.Ht,71.20.Hk 
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To improve the efficiency of polymer liglit-emitting diodes, LEDs, it is essential to un- 
derstand and improve performance of the contacts. Contact injection into LEDs made of 
MEH-PPV [poly(2-methoxy,5-(2'-ethyl-hexoxy)-l,4-phenylene-vinylene] has been attributed 
to tunneling of the electrons and holes into a broad conduction or valence band, respectively, 
in the polymer through interface barriers arising from the band offset between the polymer 
and the metal electrodes. ^'^ Good agreement of the I-V characteristic with the field depen- 
dence given by Fowler-Nordheim tunneling has been shown in references 1 and 2, although 
others did not find the Fowler-Nordheim field dependence for their diodes. 

The picture of carriers tunneling into a wide band cannot be correct for the conducting 
polymer samples now available. First, because of the short conjugation lengths for typical 
polymers such as PPV, on average ~ 6 or 7 monomers, continuum-like bands become sets 
of discrete levels. For the average conjugation length the level spacing is ~ several kT at 
room temperature.^'^ A more serious objection, however, is that injection of electrons or 
holes is into polaron levels,^'® which for electrons lie below the LUMO and for holes above 
the HOMO. The separation of a polaron level from the LUMO or HOMO depends on the 
conjugation length of the segment. For PPV calculations give the separation as 0.15 eV ^ or 
0.2eV ^ for very long segments, increasing to ~ 0.7 eV for a 3 monomer long segment. It is 
also a consequence of the short conjugation lengths that theories based on the formation of 
a bipolaron lattice in the neighborhood of the contact cannot apply to currently available 
conducting polymer samples.a 

Based on the above considerations the scenario for conduction in an LED begins with 
a carrier from the metal tunneling into a polaron level close to the contact. The carrier 
then diffuses and hops in the field from one conjugation length to another. The variation in 
conjugation lengths, and the presence of defects, result in a spread in energy of the hopping 
sites (diagonal disorder). An appropriate model for treating the transport in this system for 
low injection is the disorder model pioneered by Bassler and associates.^ The distribution 
in energy of the polaron states is taken as a Gaussian with variance aB From the expected 
spread in conjugation lengths it is reasonable that a ~ 0.1 eV.li^ 



Critical for the behavior of a contact is the location in energy of the polaron states of 
the polymer relative to the Fermi energy Ep of the metal. For specificity we will discuss 
the case of electron injection into the polymer but the results apply to hole injection with 
the usual modifications. Internal photoemission measurements, such as those of Campbell 
et al,^ yielding the energy required to inject an electron from a metal into a polymer, give 
the energy difference between Ep and some average state in the polaron distribution. We 
denote this energy by W. In what follows we calculate I-V characteristics for a case of large 
W, which means small injection, using the results of a Monte Carlo simulation based on the 
disorder modeLGJ We then carry out a calculation for W ^ 0, which is the case for calcium 
contacts on MEH-PPV, ^ using the classical approach of Rosell3 and Lampert.Ej In the latter 
case agreement is obtained, for reasonable values of the parameters involved, with the I-V 
characteristic of samples with only a Ca contact injecting, providing we take into account 
the strong variation of mobility with field documented for PPV by Karg et alM 

In the Monte Carlo simulationliil the energy f/ of a polaron site as a function of the 
perpendicular distance x from the metal-polymer interface is written 

U{x) =W - eEx - e^AKx (1) 

Here E is the electric field intensity and the last term represents the image force, k being 
the dielectric constant. In the presence of energetic disorder U{x) gives the value of mean 
energy e{x) of polaron sites at a distance x from the interface. Note that all energies 
are measured relative to Ep. With W large, e.g. 0.6 eV or greater, injection is small 
and space charge may be neglected. To calculate the incoming flux of carriers we have 
assumed that they tunnel into a polaron level with energy e at a distance x from the 
interface at the rate Vm{x)exp{—e/kT), where Vm{x) is a distance-dependent prefactor. 
With this assumption the total number of carriers tunneling into polaron levels per second 
is Vm{x) exp[{—e{x) +a'^/2kT)/kT]. Thus the energy distribution of initially populated sites 
is displaced by a'^/kT from the available site distribution. By Monte Carlo simulation we 
followed the hopping of the carriers through a sample of 12 layers, 331 x 331 sites/layer. 



finally obtaining the yield, i.e., the fraction of injected carriers that escape the return to 
the electrode and reach the opposite boundary of the samplcEJ Even in the absence of 
disorder the image force results in the great majority of carriers returning to the electrode 
at low fields. An analytic solution for the ordered case based on Eq. (1) (neglecting space 
charge )li3yj gives the yield as 0.3% in a field of 1.25 x 10^ V/cm at 300 K for the parameters 
of Ref. |ll]. A field ~ 10 times as large is required for essentially complete collection. The 
results of the analytic solution were in good agreement with the Monte Carlo simulation for 
cr = 0, indicating that the simulation sample was thick enough.EJ In the disordered case the 
carriers have to overcome the random barriers due to disorder as well as the image force. 
The result is that at fields of ~ 10^ V/cm mild disorder {a = 0.08 eV) results in a yield 
smaller by a couple of orders of magnitude even though the injection is larger. In addition 
the yield increases more strongly with increasing field.EJ Current vs. field may be obtained 
by multiplying the yield by the number of carriers entering the polymer per second. As 
shown in Fig. 1, for W = 0.7 eV and a sample length of 120 nm, the effects of the image 
force and moderate disorder are a current increase by a factor of 10^ as the applied voltage 
goes from 1 to 20 volts. 

The photoemission data of Ref. 5 give the Fermi level for Ca lying above the polaron 
level in MEH-PPV by 0.05 eV. This suggests that Ca provides an ohmic contact on MEH- 
PPV, meaning that it could supply the maximum required current, i.e., space-charge limited 



current. Calculations of current vs voltage for ohmic contacts were carried out in Refs. 12 



and |13[ The calculations are based on Poisson's equation: 

{K/e){dE /dx) = {n — fi) + {rit — fit) , (2) 

where n and rtt are the densities of free and trapped electrons, respectively, and n and fit 
their respective average values for the sample in thermal and electrical equilibrium with the 
contact (no apphed voltage). The equation for current density J is simplified by neglecting 
the diffusion terms as was done in Ref. |T^. In this case we include the field dependence of 
mobility /i found for hopping in many disordered systems, including holes in PPV,li3 to give 



for the steady current 

J = nei^oe^'^E (3) 

where Hq is the zero-field mobihty and a was taken as a parameter in the calculations. With 
the simplification of neglecting diffusion current the boundary condition at the cathode 
interface is i? = at x = 0. The other electrode is taken as non-injecting. 

Although there is an estimate of the trap density in MEH-PPV, specifically a few times 
10^^/cm'^,li3 we lack information about their location. We therefore carried out numerical 
integration of Eqs. (2) and (3) for trap-free {rit = fit = 0) and all-traps-filled {rit = Nf, the 
total trap density ) cases. As discussed by Lampert, the trap-filled case does not give the 
correct current at low fields. Instead there is a voltage threshold for current fiow because 
before voltage is applied there is already unneutralized charge in the traps which prevents 
the injection of additional charge at the electrodes. When Nt ^ n, as is likely to be the 
case here, the current rises very steeply with voltage beyond the threshold. Nevertheless the 
trap-filled solution should be good at high enough fields. 

In Fig. 2 we compare our calculated results with the experimental data of Parker for a 
Ca contact on MEH-PPV in the electrons-only case, i.e., with the work function of the other 
electrode too low to contribute significant hole current (at least below 18 V)^. It is seen 
that good fits can be obtained for both the trap-free and trap-filled calculations, although, 
as anticipated, there is no fit at the lower fields in the latter case. For the trap-free case the 
parameters for the fit shown are a = 8 x 10^'^cm-'^/^ /V^^"^ and /iq = 5 x lO^"*^^ cm^/Vs. For 
the trap-filled cases both solid and dashed lines correspond to /io = 5 x 10~^ cm^/Vs. For 
the solid line the other parameters are a = 4 x 10^^ cm^/^/V^^^ and Nt — nt — n = 6 x 10^^ 
/cm'^, while for the dashed line a = 4.5 x 10~^ cm^/^/V^/^ and Nt — fit ~ n = 10^^ /cm^. 
Smaller /iq is required for the fit to the trap free case because traps are not present to keep 
the current down at a given field. The parameter values are reasonable, particularly those 
for the solid line in the trap-filled case. Extrapolation of the fi vs E data for holes of Ref. 



T4| to -E = yields /iq = 5 x 10 ^ cm^/Vs. Electron mobility is thought to be considerably 



lower than hole mobility in PPV, the difference most likely being due to deep traps (perhaps 
carbonyls) for the electrons. It is not unreasonable that with Ca contacts the deep traps 
are filled and electron mobility becomes comparable to hole mobility. Karg et al. obtain 
a = 6 X 10^'^ cm^/^/V^/^, close to the a values obtained here. Finally, since n is expected to 
be small, Nf — fit is close to the trap density estimated by Campbell et a/ from the magnitude 
of the initial increase in capacitance with forward bias.li3 From the good fit at high voltages 
we conclude that the current there is space-charge limited current. It does not vary as V^ 
because of the strong field dependence of /i. We note that the importance of space charge 
effects has also been stressed, albeit within the Fowler- Nor dheim model, by Davids et alx^ 

In summary, we have pointed out that Fowler- Nor dheim tunneling can not describe 
contact injection into currently available polymer samples because the short conjugation 
lengths mean they cannot have broad bands. Injection and transport involve only polaron 
levels. These have a spread in energy due to the conjugation length variations and other 
defects. This disorder, even though relatively mild, can decrease injection by a couple of 
orders of magnitude and makes mobility highly field dependent. Ca contacts to MEH-PPV 
should be ohmic because Ep lies above the polaron level. I-V characteristics resulting 
from injection at Ca contacts are well fitted by theory for space-charge limited current with 
reasonable values for the mobility and its variation with field, and the trap density. 
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FIGURES 
FIG. 1. Current density vs voltage from Monte Carlo simulation for a disordered sample 

{a = 0.08 eV) with one injecting contact with large W (0.7 eV). Sample length is 120 nm. 

FIG. 2. The dots represent I-V data from Ref. 1 for 120 nm long MEH-PPV LEDs with one 
Ca contact and the other contact Mg (•) or Nd (o). The lines represent theoretical fits to the data 
for the trap-free and trap-filled cases. 



